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THE STABILITY OF ORTHOTROPIC ELLIPTIC CYLINDERS 

IN PURE BENDING-* 
By 0. S. Heck 



The theoretical critical "bending stress of elliptic 
cylindrical shells is determined on the assumption of in- 
finite shell length and absence of local instability phe- 
nomena. The results of tests on isotropic elliptic cylin- 
drical shells stressed in "bending are compared with the 
theoretical results. The practical applicability of the 
theory is discussed. 

I. INTRODUCTION 



The preliminary calculation of the load capacity of a 
thin-walled cylindrical shell under "bending stress is of 
importance in airplane statics in the analysis of shell 
"bodies. Hereby it does not merely pertain to isotropic 
circular cylindrical shells which, in the literature up to 
now, are almost exclusively used in bending tests, but also, 
above all, to the bending of orthotropic (orthogonally ani- 
sotropic** and stiffened cylindrical shells. The study of 
shells with other than circular sections (elliptic, for ex- 
ample) is of particular practical importance. 

General application of Navier's simple bending theory 
to thin-walled beams is no longer permissible, according 
to Prandtl, because of the occasionally enormous strain of 
the section under load. On the strength of this, and 
prompted by experiments on elastic balance tubes, von 
Karman investigated the bending of curved thin-walled 
pipes (reference l). But his findings are inapplicable to 
straight pipes. Besides, von Karman failed to mention the 
important fact that the bending moment has a maximum (crit- 



" ii 
* ,f Uber die Stabilitat or thotroper ellipt i scher Zylinder- 

schalen bei reiner Bie^ung." Luf tf ahrtf or schung , vol. 
14, no. 3, March 20, 1937, pp. 137-147. 
**The concept of orthogonally anisotropic plate was proba- 
bly first given by M. T. Euber (cf. Bauing. , vol. 4, 
1923, p. 354. 



2 



N.A.CA. Technical Memorandum No, 834 



ical) value conditioned "by the flattening of the cross 
section. 

The first study on the pure "bending of infinitely 
long isotropic circular cylindrical shells, was made "by 
Brazier (reference 2). He shows that the bending moment 
does not increase linearly with the curvature of the shell 
axis and that it assumes a maximum (critical) value for a 
certain curvature of the shell axis; for greater curvature 
the system ceases to be stable. 

This type of instability (similar instability pre- 
vails under compression and bending of a bar) is essen- 
tially different from that encountered in the usual sta- 
bility problems, such as in the buckling of a straight 
bar, for example. 

Figure 1 shows the behavior of a straight compres- 
sion member* The load P runs linear to the approach f 
of the bar ends as far as the branching-of f point (buck- 
ling load). From that point on, there are two possible 
conditions of equilibrium: the bar may remain straight - 
that is, in unstable equilibrium, or deflect sideways un- 
der stable equilibrium. In both cases the load P in- 
creases beyond the buckling load, 

Erazier f s discussion of bending of a cylindrical shell 
is illustrated in figure 2, The curve which gives the 
bending moment B relative to the curvature K of the 
shell axis, has no branching-of f point. The bending moment 
cannot be increased beyond the critical value B^ r , After 
the critical condition is exceeded there remains only one 
possible equilibrium condition of the shell, and that is: 
unstable. 

Figure 3 finally shows the behavior of a cylindrical 
shell in bending when accompanied by local buckles which 
belong to the instability phenomena of the usual type. 
The maximum load supported by the shell is more or less 
reduced by the buckles; hence Brazier's value obtained by 
disregarding a local buckling, constitutes a theoretical 
upper limit of the maximum load in an infinitely long 
shell . 

II, FORMULATION OF PROBLEMS AIJD ASSUMPTIONS 



In the present report the behavior of orthotropic 
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cylindrical shells of elliptic section in pure "bending is 
theoretically investigated. Critical "bending moment, 
critical bending stress, a-nd deformation of shell section 
are determined. The principal assumptions are as follows: 

1. The cylindrical shell is of infinite length. 

2. There are no local instability phenomena (i.e., 

no wrinkling). 

5, The stresses remain below the proportional limit 
of the material. 

The value of our findings for practical application, 
if premises 1 and 2 are not met, is discussed in the next 
section. Whether assumption 3 is fulfilled can always be 
vcrif i od» 

The locrl instability of long shells can "be investi- 
gated when allow&nce for the strain condition - i.e., the 
change in curvature radius of the shell section and of the 
distances from the zero line - is made in accordance with 
the considerations of the present report. An approximate 
study of local instabili ty^ has been made "by Brazier for 
the case of isotropic circular cylindrical shell. But the 
local instability is so seriously affected by inevitable 
preliminary wrinkling that theoretical studies which do 
not allow for these preliminary wrinkles, are of limited 
practical importance only. 



The critical bending moment B kr of an orthotropic 
cylindrical shell of elliptic section stressed in bending 
in the major or minor axis of the section is on the assump- 
tions cited in the preceding section: 



III. RESULTS AND SCOPE OF VALIDITY 



B kr ^ W CT kr 



(1) 



with 




Hereby 
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E is the modulus of elasticity of the material 

U,, Poisson's ratio (about 0,3 for steel and duralumin) 

s, wall thickness of shell 

s m> . mean wall thickness, which governs the stiffness of 

the shell to axial stresses (in the case of the shell 
with closely spaced stringers, which approximately 
resembles an orthotropic shell, it is: s m = s + 

3? 

where F = total section of stringers, and u * 
circumference of section of mean shell surface)* 

a,b, major and minor half axis of the shell section 

p, radius of curvature of undeformed shell section at 
the point of stress peak (p = a 2 /h under "bending 

on the major sectional axis = "b / a by bending on 

the minor sectional axis) 

W, section modulus of undeformed cross section of the 
elliptic cylinder 

(W g = J b (b + 3a) s m 
referred to major sectional axis 

W k = | a (a + 3b) s m 
referred to minor sectional axis) 

°kr > critical bending stress (cr kr is simply a fictitious 

stress, because W refers to the undeformed cross 
section of the shell) 

c, a numerical value which depends on axes ratio a/b 

2 a 2 - b 2 

of the shell section or the value k = 5 — and 

a 

may be obtained from figure 4. 

The subscripts g and k indicate, respectively, the 
bending about the major and minor cross-sectional axes. 



*Uniform stringer distribution being assumed. 
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For a shell having unlike moduli of elasticity in Ion 
gitudinal and circumferential directions, the critical 
"bending stress is: 



The "real 11 critical "bending stress oy_ r 1 ( in contra- 
distinction to the "fictitious" critical stress (Jkr ) must 
not - if the findings of this study are to retain their 
validity - exceed the proportional limit of the material* 
To decide whether this condition is fulfilled, compute 
a^ r 1 according to (2) or (3) whereby c g 1 and c k f ob- 
tainable from figure 5, replace Cg and c^. 

The extent of the cross-sectional strain may he seen 
from figures 6 and 7, where the relative length changes 
in the axes of the shell section, on reaching the critical 
condition, are plotted** 

In the calculation of the quantities plotted in fig- 
ures 4 to 7, the squares and products of the displacement 
components of the shell element in their cro ss-sect ional 
plane were disregarded relative to the first powers of 
these values, as in Brazier*s study. To gain an idea of 
the effect of this omission, we made a more accurate cal- 
culation for k2 = 0 and k 2 = 0.3, the results of which 
are also included in figures 4 to 7. But oven these val- 
ues are not absolutely correct, since there are still 
other influences which are not considered in the calcula- 
tion. In the calculation of the curvature change of the 
shell section, for example, the displacement components 
of the shell elements in the cross-sectional plane are 
considered as small; i.e., higher powers of these values 
are neglected relative to the first powers. Further, it 
was assumed that the strain in the shell section consists 
solely in a deflection, but not in a length change of the 
line elements of the shell section. The more exact calcu- 
lation gives greater values for o^r • s0 that ^e app^ 1 "* 



w is the shifting of a cross-sectional point at right 
angles to the circumference (positive inward). The expla- 
nation for w 0 and wrr/2 will be found in figures 12 




(3) 



where E 1 = modulus of elasticity in circumferential, and 
S 2 in longitudinal, direction of the shell. 



and 13. 
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cation of the results of the simple calculation leaves 
one on the safe side. Experimental verification of the 
theory seems, in the face of these facts, absolutely nec- 
essary. 

In the following, the extent of the practical applica- 
tion of equations (1) to (3) to cases where assumptions 
(1) and (2) are fulfilled, is discussed. 

a ) Isotropic shells,- In isotropic shells, even if 
of great length, the "beginning of the instability through 
collapse of the shell- walls toward the neutral axis, is 
initiated "by a local wrinkling on the compression half of 
the cylinder, as a result of which the failing load of the 
shell is reduced, Numerous experiments on isotropic cir- 
cular, cylindrical shells (cf. figs. 17 to 19) have shown 
that the average maximum "bending moment supported "by the 
shells is not very much different from the theoretical 
value for the infinitely long shell without consideration 
to local instability phenomena. 

The average value of c of 77 tests is c a verage = 
0.357. But the scatter of the experimental values is in 
part quite considerable. From the available test data 
there is no indication of any effect of shell length (fig. 
19). 

For the purpose of checking the applicability of the 
theoretical values arrived at with the assumptions out- 
lined in section II, for elliptic cylindrical shells of 
finite length, we made several experiments with isotropic 
elliptic cylindrical shells of varying axes ratios. The 
tests disclose through the theoretical values computed 
with the coefficients Cg and c^ a good agreement for 
bending a.bout the major cross-sectional axis, but markedly 
lower experimental values for bending about the minor 
cross-sectional axis, (Cf, section IV and fig. 4.) 

In the latter case, it means that local wrinkling of 
the shell wall results in a material reduction of failing 
load of the shell, contrary to the experiences with circu- 
lar cylindrical and elliptic cylindrical ■ shell s in bend- 
ing about the major cross-sectional axis. A certain ex- 
planation for this fact may be found when assuming that 
the failure then occurs in the event that the axial stress 
0" in the shell reaches at some point of the cross section 
the critical value ok r of equation (2). This most un- 
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favorably stressed point follows from the condition that 
the product of curvature radius of shell section and the 
distance from the neutral zone must "become a maximum. 
When disregarding the strain of the shell section this 
point is given by 

t = arc sin — - (k 2 > 0,25) 

Then the critical "bending moment can "be computed from the 
formulas (1), (2), and (3), is at k 2 > 0.25 instead of 
c^- the coefficient 



3 /2 

k (1 ~ k 2 ) c k 



^ k ~ 3 JS 

is substituted. At k 2 % 0.25 the extreme fiber on the 
compression half of the shell is most adversely stressed. 
Figure 4 shows c k plotted against k . The agreement 
with the experimental values is comparatively good. 

b) Orthotropic shells,- While the local bulging of 
the walls of a shell with lengthwise closely spaced stiff- 
eners, which approximately resembles an orthotropic shell, 
is dependent on the bending stiffness of these stiffeners, 
the failure due to collapse of the shell walls is depend- 
ent on its cross section. There is no local bulging of 
shell walls before failure in very long elliptic shells 
with closely spaced stringers for sufficiently high bend- 
ing stiffness and small cross section. In that case the 
results of the present report are exactly valid. 

Ar. concerns the applicability of the theoretical for— 
mula and the effect of shell length in orthotropic ellip- 
tic shells (plywood shells, for instance), no experiments 
pre available. For computing the failing moment of short 
stiffened shells with strong fairly closely spaced frames, 
the results of this study are inapplicable as the re- 
strained flattening of the shell section increases the 
failing load of such shells as a rule very considerably 
compared to the failing load of very long shells without 
frames.* If the strain of the section of such a shell is 



*The theoretical formula (2) is also useful for calculat- 
ing the closed part of a shell body with comparatively 
closely spaced intermediate bulkheads, since it indicates 
the critical bending stresses for the collapse of the shell 
walls in buckling form; that is, a lower limit in shells of 
(Continued at bottom of page 8) 
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negligibly small, then the study of the local stability of 
the assumedly orthotropic shell is comparatively simple as 
against the general case, "because the "bending moment up to 
stability limit runs linear with the curvature of the shell 
axist As regards the bulging of the walls of orthotropic 
shells, a better agreement between theory and test seems 
indicated than with isotropic shells, since the inevitable 
preliminary bulges are smaller in comparison with the mean 
wall thickness. 



IV. EXPERIMENTS WITH ISOTROPIC ELLIPTIC CYLINDERS 

A series of failing tests in pure bending about the 
major and r.inor cross-sectional axes was carried out on 
isotropic elliptic cylindrical shells of duralumin of two 
different axes 1 ratios. The dimensions of the test speci- 
mens are given in table I. (I is the free length of the 
cylindrical shell.) The modulus of elasticity of the ma- 
terial was established at E = 7.5 x 10 5 kgcnr 2 . The 
evaluation of the tests was made with with v = 0*3 
Poisson*s ratio. The experimental arrangement itself is 
illustrated in figure 8. The wooden frames at the ends of 
the test cylinder serve to press the sheet against them 
through exactly fitting wooden jigs. One en" frame is 
clamped to a solid frame while a pure bending moment is 
applied at the other frame. The weight of this frame is 
compensated. Two failing tests could be made on each 
specimen, After the first test the cylinder was turned 
through 180° on its axis, which left the still undamaged 
part of the specimen on the compression half of the cylin- 
der for the second test. 



TABLE I, Dimensions of Test Specimens 



Sp ecimen 


a 


■b 1 


' I 


s 


cm 


cm 


cm 


cm 


1 


22.5 


15 


~78 


0.050 


p 


SO 


15 


127 


.051 



* (Continued from footnote, page 7) 

of finite length, which are bounded by fixed ribs. Logic- 
ally, the thickness of a smooth sheet must thereby be writ- 
ten for s in (2), whose bending stiffness equals the mean 
bending stiffness of the stiffened shell in circumferential 
direction under consideration of the frames. One obtains 
in this manner a control point for the rib dimensions. 
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The results of the failing tests are given in table 
II, For comparison with the theoretical values (admitted- 
ly derived for infinitely long cylindrical shell and dis- 
regarded local instability phenomena), we have included 
the experimental values of Cg and c^ in figure 4. A 
discussion of the test data is given in the preceding 
section. 



TABLE II, Test Data 





Bending about major 
cross-sectional axis 




Bending about minor 
cross-sectional axis 


Spec- 
imen 




ff_ 

6 


s kr 


e kr 




Pk 


^k 


B k 1 
kr 


K kr 


C k 




cm 


cm 3 


cmkg 


kgem" 3 




cm 


cm 3 


cmkg 


kgem -2 




1 


33.75 


48.5 


19600 
21200 


405 
437 


. 347 
.375 


10 


59.6 


42500 
47300 


713 
794 


.181 
.202 


2 


60 


63.1 


17350 

16900 


275 
-2*8 


.411 
.401 


7.5 


90.1 


5 3700 
52200 


595 
580 


.111 

.108 



V. THEORETICAL ANALYSIS 
1, formulation of Problem as Variation Problem 



To analyze the behavior of a circular or elliptic cyl- 
indrical shell in pure bending, we resort to the principle 
of minimum potential energy. The assumption of infinite 
shell length neutralizes eventual edge effects. Conforma- 
ble to another assumption, there is to be no local bulging 
of shell walls. Then all sections of the shell are strained 
in the same fashion for reasons of symmetry, so that the 
analysis can be restricted to a piece of shell of length 1« 

For a given curvature ft of the shell axis, the form 
which the shell section assumes must be so defined that by 
fixed K the energy of form change U becomes a minimum. 
So if the strain of the shell section and consequently, 
the energy of form change U is known in relation to cur- 
vature ft, ' the bending moment 3 follows as derivative 
of U with respect to k. 
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The form change energy U comprises: 

1) The form change energy tlx % corresponding to the 

strain of the shell section. 

2) The form change energy U 2 , corresponding to the 

length changes of the shell fibers. 

The strain of the shell section is assumed to consist 
only in a deflection hut not in a length change of the lin- 
ear elements of the section of the middle of the shell. 
It is assumed that both strains corresponding to the energy 
of form changes JJ 1 and U 2 take place successively, for 
example - first the strain of the shell section and then 
the elongations or contractions in shell-length direction. 
The additional strain of the shell section following the 
initiation of the axial stresses as a result of transverse 
contraction, is small enough to be negligible. 



It is then: 



E 



J m K 2 J 2 (5) 



wi th 



J x = § U K U ) 2 du (6) 
J P = § h 2 du (?) 



'2 



whereby the integrals extend over the whole circumference 
of the shell section. It is: 

K, curvature of the shell axis 

&K V , change in curvature of the section of the median 
shell area 

u, arc length of the section of the median shell area 

h, distance of an element of the strained shell from 
the neutral axis 

Assume the strain of the shell section to be defined 
by the natural coordinates v and w of the displacement 
of the shell elements in the sectional plane: v, to denote 
the component of displacement in tangential direction (pos- 
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itive in the sense of increasing u) and w, the component 
in normal direction (the inside normal to "be positive) of 
the undeformed section. The curvature change A K u of 
the shell section, which may he any, so long as it has no 
corners, can he expressed with Frenet 1 s formulas through 
the displacement components v, w, and their derivatives* 
We have: 

w P 1 , V 

A K n = -o- + w n ~ v -o (8 ) 

P P 

Here p is the radius of curvature of the shell section 
(positive if the center of curvature lies on the inside 
normal), p ! the first derivative of p with respect to 
arc length u, and w" the second derivative of w with 
respect to u. If, as according to our premise, the length 
of the linear elements of the section is constant, then 
there exists "between v and w, the relation: 

w = p v» (9) 

With due regard to (9), equation (8) "becomes: 

A k u - p v m + 2p « v" + (p n + v f - v (10) 

It will he ohserved that in (8) to (10) the higher 
ipov^er^ of v and w and their derivatives are disregard- 
ed relative to the first powers of these quantities, as a 
result of which the validity of these formulas is confined 
to relatively small v and w. 

The prohlem now is to so define v and w that for 
given curvature K of the shell axis the energy of form 
change U assumes a minimum value; i.e., to solve the 
variation prohlem: 

U = U 1 + U 2 = Min (11) 

2. The Orthotropic Circular Cylinder Under Pure Bending 

We first analyze the pure "bending of the infinitely 
long orthotropic circular cylinder on the premises of 
section II, and solve the previously derived variation 
prohlem for this particular case. The curvature change 
Aku 01 ~ the shell section "becomes: 



A k u = r v ,! ' + 2- 

r 



(12) 
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according to equation (10), with r equal to radius of 
median shell area. The distance h of one element of the 
strained shell from the neutral axis is: 

h = r cos t - v sin t - w cos t (13) 

wherehy w is to "be replaced according to (9); the sig- 
nificance of t is seen from figure 12. The variation 
problem (11), for which the solution could equally well "be 
arrived "by integration of the correlated Eulerian differ- 
ential equation, is solved directly "by the formula (Ritz's 
method ) : 

n 

v = r Z Ai sin 2 J t (14) 

IT o terms other than those given can appear in the formula 
for v for reasons of symmetry. With (14) and allowance 
for (12), (13), and (9), equations (6) and (7) give: 

J, = f £ f (4 I* - I)' A; 2 (15) 
J 2 = tt r 3 (1 - 3 A x ) (16) 



In the determination of J 2 the squares and prod- 
ucts of v and w are neglected relative to the first 
powers of these values. The conditions for minimum energy 
of form change U are the n equations: 



or 



with 



3A 



!i- = 0 



a j 



3A 



( j = 1 , . . . , n ) 



l _ 



12 n a_j £ 

r 4 8 A"; 



s m k 3 r 



3 



(17) 
(18) 



With (15) equation (17) "becomes: 

9 J" (M « - if Aj . - || 

From these equations, the quantities Aj can "be computed: 
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Hereby 11 may "be of arbitrary size; that is, 

v = r s i n 2t 

2 

is the exac-t solution of tire variation problem (11). The 
"bending moment B "becomes: 

B s M = |S s E s m K J a = tt E s m k r 3 (l - ^ (19) 
dK 3k m 2 \ 2 / 

The "bending moment reaches its maximum value B kr if 

equation 

11 = d B + 8 B dN = 0 (20) 
dK 3 k 9B dK 

2 

is fulfilled. Equation ( 20 ) gives: II = — . 
By observing (18), we have: 

B kr = C --g r s jnr n (21) 



u 2 



0 = = 0.987 

9 

Putting B kr = W C7 kr 

whereby W = tt r r 

we find: 

E s 

a - c — 

r- r r a/ s T 



5 m 



(22) 



2 v ( 



2.V.2 = 0.314 



but a kr is only a fictitious stress since T7 is the sec- 
tion modulus of the unstrained circular section. 

Allowing for the squares and products of v and w, 
disregarded in equation (16), we have: 
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2 ' 

J 2 = rr r 3 (1 - 3A X ) + jf (v sin t + w cos t ) ds 



= tt r» (l - 3A X -1- | A x s 



A, A : 



17 



45 



37 



• ~ A 3 " + 3* A 2 A 3 + ^ A 3 C 



...) (23) 

Then the equations (17) read, for n = 3, for example: 
3 Ai + 



§ Ai + ? Aa^) N 
2 4 / 



300 A s + 



3675 A* + 



(I 



Ai + ^ A 2 + 



¥ A 0 



ft 5 ^ - f a,) 



N = 0 



N = 0 



(24) 



from which Ai to A3 (up to A n in general) can "be ob- 
tained. The tending moment 



B 



dU _ oU _ 



- = E s m K J 



dK dK 

reaches a maximum value if 
dB 

dK dK ' dK j=i dAj dll 



dB dN £ 9b dAj _ 



(25) 



(26) 



This equation gives H and, "by observing (18), (23), 
and (25), the critical "bending moment. It again yields 
the equations (21) and (22) for 3 kr and a^ r hut with 
different coefficients C and c. It is: 

C a 1 # 22 and c = 0.388 

which values are, as stated before, still not completely 
correct • 

On reaching the critical condition (fig. 12) the rel- 
ative shortening in diameter of the shell section perpen- 
dicular to the neutral axis is: 



— — = N = - 
r 9 



(27) 



if the squares and products of v and w are neglected 
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in (7). The relative lengthening of the diameter coinci- 
dent with the neutral axis is of the same magnitude. If 
the squares and products of V and w in J 2 are con- 
sidered, then the relative shortening of the diameter per- 
pendicular to the neutral axis amounts to 



* = 2 X i Aj (28) 

and the relative lengthening of the diameter coincident 
with the neutr a 1 ax i s 

- Mi = 2 § + 1 J Aj (29) 

r j=l J 

Hero the values A,- should be determined from (17) and U 

fron (26). The values w 0 /r and - % / g/ r together with 

the correlated values of shells with elliptic section are 
shown in figures 6 and 7. 

The actual critical "bending stress cr^r T contrast 
to the fictitious "bonding stress CTw) "becomes: 

^kr' = fT 

Thereby the critical "bending moment is, according to (19): 

B kr "-- E s m K J a 

and the section modulus of the deformed shell section (re- 
ferred to the neutral axis): 

W* = — 

Observance of (18) gives: 



r - w 0 



E s / s 

°kr " 



| = C ! 



J 1 - V 



2 r 



If the squares and products of v and w in equation 
(7) arc neglected, it affords with equation (27): 

c 1 = ( 1 - IT ) f$ 
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Allowance for the squares and products of v and w, and 
observance of (28) gives: 

o' « <1 - 2 jit J Aj) JT 

The values computed for c 1 are plotted in figure 5. 

3, Pure Bending of the Orthotropic Elliptic Cylinder 

In the following we solve the variation problem (11) 
for a cylindrical shell of elliptic section. The study is 
confined to the symmetrical cases of bending about the ma- 
jor or minor axis of the shell section. A rectangular sys- 
tem of coordinates in the plane of the shell section makes 
the x-axis coincident with the major, and the y-axis coin- 
cident with the minor, cross-sectional axis (fig. 13). 
Then the equation of the section of the median shell area - 
expressed in parameters- reads: 



x = a sin t 



y = b cos t 



where a is the major, and b the minor, half axis of 
the cross-sectional ellipse. Conformable to (6) and (10), 
we have: 

2 



whereby 



+ 2p 



! 



! 



p" + -)V« - V -3 



du 



(30) 



du 



a (1 - k 2 sin 3 t) 1/ ' dt , 



, s a 2 - "b 2 
k = 



S (1 - k 2 sin 2 t) 3/2 a, 

D 



p ' = - | | k 2 sin 2t 



p» = - 3 | k 2 cos 2t (1 - k 2 sin 2 t)~ 1/S \ 

1 o .a 

If h g is the distance of a shell element from the 
neutral axis in bending about the major cross-sectional 
axis, and h k , the corresponding distance in bending about 
the minor axis, we have (cf. fig. 13): 



h g = b cos t - w cos cp - v sin cp 



h 



k 



sin t - w sin cp + v cos cp 



(31) 
(32) 
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Here cp is the angle of the normal with the minor axis 
of the ellipse. It is: 

sin cp = £ sin t (1 - k 2 sin 2 t )~ 1 /2 
a 

cos cp = cos t (1 - k 2 sin 2 t)"* 1/2 

According to equation (7) the observance of <3l) and 
(32) in "bending about the major cross-sectional axis (de- 
noted hereafter by subscript g) gives: 

J 2 = / ("b* cos 2 t-2b w cos t cos cp -2b v cos t sin cp)du 

5 (33) 

and by bending about the minor cross-sectional axis 
( subscript k) : 

J 2 = / (a 2 sin 2 t-2a w sin t sincp+2a v sin t cos cp)du 
k (34) 

The squares and products of v and w are neglected 
against the first powers of these values; p v 1 can be 
substituted for w according to (9). 

The variation problem (11) can be solved by the Ritz 
method. The function v can be approximated by 

n 



v a a .Z Aa sin 2j t (35) 



For reasons of symmetry, no terms other than those 
given can occur in the formula for v. With formula (35), 
equation (30) becomes, after several intermediate compu- 
tations: 



tt/2 



Jl = a (1 - k 2 ) f 



(jla A J K j) dt (36) 



wi th. 



= |2j cos 2j t [-4j 2 +(a+p sin 2 t+Y sin 4 t) 

(1-k 2 sin 2 t)" 2 ] + e sin 2t sin 2j t(l~k 2 sin 2 t ) _1 



[4j 2 + £ (1-k 2 sin 2 t) S ]| (1-k 2 sin 2 t) 1 
where for abbreviation: 



/4 
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1 - 3k 2 = a, 4k 2 = 8, -2k 4 = Y 

| k 2 = € , 1 - k 2 = t 

The term for Ji given in (36) is valid in lending 
ahout the major, as about the minor, axis of the shell 
section. The evaluation of the integral in (36) "by series 
development "being very tedious, the numerical integration 
is carried out "by the G a us s-Lobat t o method (reference 8). 
The integral is approximately replaced by a sum of n 
terns. The summands are the values multiplied "by certain 
weights g of the to-be integrated function at the inter- 
val stops and at m - 2 prescribed points of the interval 
The thus-obtained approximation is of the order of 2m - 1 
i.e., a parabola of the 2m - 1 degree is exactly inte- 
grated hy this method. Then equation (36) becomes: 

j a £IL_ E g ( r A, Kj J (37) 

1 a (1 - k 2 ) p = i P VJ«i J J P' 

The integrals J P and J a , , conformable to equations 

g * 



J 2g 



(33) and (34) can he exactly defined, i.e.: 

= 2tt a 3 j^(l-k 2 )R--' 3 -™^ A, - A 

[ 

(l-k 2 sin 2 t) 1 2 cos 2 t at = I - \ (I) k 
6 \2«4/ 3 8 V2*4«6 



(38) 
(39) 



■with. 



R = — 

TT 



/ 5 

S = | ^(l-k 2 sin 2 t) 1/2 sin 2 ttt-J.J (I)" k 8 

0 5 (±-_Z\ 1st „ 1 (liSlSt k! - ... 

6 \2»4/ 3 8 \2-4-6/ 5 

With (35) the energy of form change U "becomes a 
function of the n variables Aj (j = 1, B ) . and 

so variation problem (11) reduces to a common "extremum" 
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problem, The accuracy of the method can he raised at will 
"by increasing n. But, in general, it affords no possibil- 
ity to assess the errors. The values A ^ , for which U 
becomes a minimum by constant K , follow from the n lin- 
ear equations: 



o r 



£2- = o 



( j = 1 , .... n ) 



12 N 3J 2 
a* dAT 



with 



s m K 2 a 4 (1 - v 2 ) 

n = -a ^ 



(40) 
(41) 



From equation (37) follows 



d Ji 



4tt 



m 



SAj a (1 - k 2 ) 



2 Ai Z gt) (E, K ± ) 

1 = 1 TVSl P ^ 1 



Then the equations (40 ) with observance of (38) and 
(39) become, in bending about the major cross-sectional 
axis : 



.£ Ai 2 so < K i K i)p = 3 (1 - k*) (3 - k") N 

1 = 1 p = 1 3f & 

. S Ai S fo (Za Ki) p = 3k 2 ( 1 - k 2 ) 8 

1=1 1} " 1 1 - 



1 = 

and in bending about the minor cross-sectional axis; 



(42) 



. S Ai Z g p (K x K ± ) v = -3 (1- k 2 ) (3- 2k 2 ) N 

1 = 1 pal 

1 Ai g g v (Ka Ki)^ = -3k 2 (1 - k 2 ) I 

1 Ai 1 g v (Kj gf)^ =0' <j - 3, ...,n) 

l — i T) — 1 



(43) 



Having defined the values Aj from (42) and (43), the bend- 
ing moment 3 follows from 
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dK 3k m 2 



(44) 



with J 2 as given in (38) and (39). The bending moment 
B reaches its highest (cri%ical) value Bj^ , if 



dB 8B 
dK ~ 3k 



+ «Y2L + 2L o (45) 

dK vdAi dN 8A 3 dN / 



From these equations follow: 



IT, = 



_2 (1 - k a ) R 
(3 - k 2 ) |± + k £ 



- 2S 



N 



(3 „ ok 2 } — + k 2 

whereby in the expression for N the values A x and A 2 

are to be taken from (42) and in the term for N k from 
(43). 

By observing (38), (39), and (41), equation (44) gives 



with 



B kr - C 



E 



— as/ 



s s 



m 



(46) 



4 

S 



TT (1 - k") R f 



<3 - k 8 ) 4 + k 2 S 



4 

— IT S 

e 



- 6S 



7 (3 - 2k 2 ) — + k 



a As 



with 
whereby- 



Defining a fictitious critical bending stress a kr 

3 kr = W CT kr 

W k = j a (a + 3b) s m 



'm 



N.A.C.A. Technical Memorandum No. 834 



21 



are the section moduli of the shell section with respect 
to the major and minor cross-sectional axes, we find: 



with 



c ~ = 




a kr = c g £ /JL (47) 



16 R / 6 (1 - k 2 ) R 



Ci 



16 (1 - k 3 ) S / 



9 1 + *Jl Z k 2 / (3 _ 2k 2 ) |* + k 2 Tr 

Figures 14 to 16 show the values Cg and C k and 

1 p 
™ slotted against k rf and n. A tjroof of the converg- 
c k 

ence of the method is withheld in the present report. 
Even so, the diagrams manifest the good convergence for 
small k 2 , while for higher k 2 , it is less good. But 

even in this range the curves toward which the values C 

1 . . s 

and — strive, can he plotted with sufficient accuracy. 

Ck , 
Figure 4 shows the values Cg and c k plotted against k 

Observance of the squares of products of v and w, 
disregarded in (3?) and (34), reveals the right-hand sides 
of equations (38) and (39) augmented "by the terms: 



/ (w cos qp + v sin cp ) 2 du = 7-^7-2 / ( L j ) cLt (48) 

and 

§ (w sin cp ~ v cos cp f du = a 3 / ( t Aj Ua ) dt (49 ) 



Her ehy 




= 2j (l~k 2 sin 2 t) 3/4 cos t cos 2j t 



■—1/4 

+ (l-k B ) (1-k 2 sin 2 t) sin t sin 2j t 
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2 2 3 / 4 

Mj = 2j (l-k sin t) sin t cos 2j t 



- (l-k 2 sin 2 t) " cos t sin 2j t 

Evaluating the integrals in (48) and (49), conforma- 
ble to G-aus s-Lobatto , the equations (38) and (39) are re- 
placed by 



-1/4 



= 2tt a" 



(l-k*) R - Ai - ~ A 2 



1 m / n \ d 1 



(50) 



J P = 2tt a" 



v 2 



m > n N a -I 

+ f i, % Q|, a j %) J < 51 > 

In "bending about the major cross-sectional axis, the 



equations defining Aj are: 



n 



m 



+ 12 H . £ Ai S g~ (L x L4) =3(l-k ft ) (3~k*) N 
j. = i p = 1 " P 



m 



^ (52) 



12 N Z An Z g^ (L 2 Li)^ = 3k (l-k ) IT 
i = i ~ r> = i n P 



n 



m 



1=1 1 p«i ^P v J * 'p 



m 



+ 1?. N 4 2, Ai _2 (LjLi) p - 0 (j=3,...,n) 



1 = 1 



> = i 



and for tending about the minor cross-sectional axis: 
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n m 

2 Ai Z g v (K x Ki) 
1=1 p=i 1 

n 

i-2 



+ 12 (l~k 2 ) U p Ad 2 (MiMi) 
1 = 1 p = i - * 

= - 3 (1-k 2 ) (3-2k 2 ) N 

n in 

.£ Ai E g p (K 2 Ki) 
1 = 1 p = i - L 

n m 

+ 12 (1-k 2 ) N £ Ai 2 £p (Ma Mi )-p 
i = i p = i r * 



(53) 



= - 3k 2 (1-k 2 ) N 



n 



m 



£ Ai E g p (Kj Ki) p 

1=1 0=1* 



n m 

+ 12 (l-k S ) H £ Ai E (M*Mi ) = 0 

i = i p = i * J 



(j=3 , ... ,n ) 



The "bending moment B follows from equation (44), 
whereby J 2 is given in (50) and (51). The value of S 
or IT, for which B reaches the critical value B^r > is 
computed from the equation: 

dB 3B dB g d3_ = 
dK ~ 3* d* 0 = 1 3Aj dB 



or 



2 (1-k ) R- (3-k ) Ai - k A 2 + 



1-k 



m / n \ 

E g_ ( £ An-Li ) 



- 2N 



2N dAi i2 dA 2 
> dT + k dF" 



8K n n m &A^ 



n 

for bending about the major axis and 



(54) 
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2S + (3~2k ) Ai + k 2 A 2 + 2 Z g^ [ Z JLa Ua ) 

p = i P \j=i J J P/ 



+ 2N 



2n dJLi o dA 2 
(3-2k 8 ) ^ + k s ^ 



n n m dA, 

+ 81T ill ill A J p£» g P (MiM J } P IlF = 0 < 55) 

for "bending about the minor axis of the shell section. 

The results for and 0"kr are again the equa- 

tions (46) and (47) hut with different coefficients C 
and c« The thus-obtained values Cg and c^ are in- 
cluded in figure 4. 

The strain in the originally elliptic section is ob- 
tained from the functions v and w known after the val- 
ues JLa have been determined. Assume that the relative 
length changes of the half axes of the shell section upon 
reaching the critical conditions are given (fig. 13). In 
bending about the major axis of the section, the relative 
contraction of the minor half axis b becomes: 

w P n 

~ ■ ---it .2 3 Aj (56) 
and the relative lengthening of the major axis a 



a J — 1 



(57) 



In bending about the minor axis of the section the 
relative shortening of the major half axis a becomes: 

= - a/TTI? E (-l) j + 1 j Jla (58) 
a j=i 0 

and the relative lengthening of the minor half axis b is 

« n n 



b 1 - k j = i J 



(59) 
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If the squares and products of the displacement com- 
ponents v and w are neglected against the first powers 
of these values, it is necessary to write: 



N 



g 



2 (1 - k ) R 

In y. s \ Ai .0 As 

(o - k ) - r + k y 



- 2S 



N 



N 



in the equations (56) to (59), while if these squares and 
products of v and w are allowed for, Ng and lT^ must 
be determined from (54) and (55). The numerical values of 
w 0 /h, etc., are illustrated in figures 6 and 7 for vari- 

■ 2 

ous k . 

The actual critical "bending stress CTfc r 1 (in contrast 

to the fictitious stress a^ r referred to the undeformed 

section), is readily obtainable. In bending about the ma- 
jor axis of the section, it is: 

a ' - 
g kr i 

Hereby the critical moment is, according to equation (44): 

B _ = E o_. K J" p 

g^r m ^ g 

with 

s 



K = 



a 2 yr 



- tf a v s m 



N 



g 



and the section modulus of the strained shell section (re- 
ferred to the major axis) becomes: 



<v g 



b - W n 



With observance of equation (56), we find: 
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'g 



E _s_ 

7rr^ P g 



2 n 

l g- 2 

1 - k j=i 



s 



Accordingly, the actual critical bending stress in bending 
about the minor cross-sectional axis is established at: 



E 



*J 1 - v z p k 




r 



n 



c k ' = (1-k 2 ) I 1 + 2/l - k s i x (-l)j + 1 3 Aj 

The values of Cg ! and 1 are shown plotted against k 

in figure 5. 
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VI, APPENDIX 



Results of Previous Studies 

a) Theoretical.- Brazier (reference 2) treats the pure 
tending of a circular cylinder of infinite length, wherein 
he neglects the higher powers of the displacement compo- 
nents v and w relative to the first powers of these 
quantities, 'He arrives at a differ en tial e quation which 
must he exactly solved, Chwalla ( r ef er ence 4) does not 
make this omission in his analysis of the work of form 
change corresponding to the tension stresses. His solu- 
tion is an approximation evolved on the "basis of an ellip- 
tic formula for the shape of the strained section of the 
circular cylinder. Nothing can he said about the magni- 
tude of the errors. Chwalla obtains a somewhat different 
numerical factor in the formula for the critical bending 
moment from that given in the present article - prohahly 
due to the fact that he does not use as mathematical ex- 
pression for the assumption of a constant arc length of 
the shell section the differential equation (9), which is 
only valid for small v and w. He rather defines corre- 
lated pairs of diameters of the elliptically deformed sec- 
tion, so that the circumference of the cross section re- 
mains constant and equal to the circumference of the medi- 
an surface of the undeformed shell, regardless of the mag- 
nitude of the strain. 

The principal results of Brazier and Chwalla have been 
tabulated in table III. 

Brazier likewise approximated the critical bending 
moments at which local instability phenomena (bulging) are 
imminent. Because he assumes the whole shell to be as ad- 
versely stressed as the extreme fiber in the compression 
zone, his values for the theoretical critical moment are 
too low. For this reason, and in consequence of the great 
influence of preliminary wrinkling through- whi ch the mo- 
ment, at which a perceptible wrinkling actually starts is, 
under certain circumstances, markedly reduced, the practi- 
cal value of this analysis is less great. Besides, he 
uses a formula for the critical stress of an axially com- 
pressed cylinder containing the factor r-T~^ (Southwell , s 

method), which is now omitted after the works of H. Lorenz 
and Km v. Sanden (reference 9). Without this factor, the 
valf.cs for the critical bending moment become slightly 
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greater. But the noteworthy fact is, that this calculation 
of the critical "bending moment results in a formula of the 
same construction as Brazier's other consideration. 

h) Exp o r i m en t a 1__ da t a . - The available results of tests 
on circular cylinders in "bending are plotted against r/ s 
and l/r (I « length of cylinder) in figures 17 to 19 
(references 2, 5, 6, and 7). 



TABLE III. Theoretical Results "by Brazier and Chwalla 





According to 
Brazier 


According to 
Chwalla 


1. 


Critical tending moment 


E 

B kr =0.987 ; o r s 2 








(with u=0.3) 


B kr =I # 035 E r s 2 


B kr =1.19 E r 


s 2 


2. 


Pertinent curvature of 
shell axis (v=0.3) 


K = 0.494 r S 2 


K=0.806 -^2 




3. 


Per tin en t relative 

shortening of diame- 
ter perpendicular to 
the neutral axis 


= 0.222 

r 


-f = 0.365 




4. 


Pertinent relative 

lengthening of diam- 
eter coinciding with 
the neutral axis 


= o.222 

r 


-H/i si 0.307 
r 




5. 


Fictitious critic al 
"bending stress 
(u=0.3) 


(T kr = 0.329 E f 


0" kr = 0.379 E 


r 


6. 


Real critical "bending 
stress (u=0.3) 


a kr ' = 0.385 E ~ 


a kr T = 0.523 E 


s 
r 



Translation "by J. Vanier, 
National Advisory Committee 
for Aeronautics. 
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Fig** 1,2,3,4,5,6 
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Figures 1-3.- Instability phenom- 
ena on elastic bodies. 
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Figure 4.- Values c' f c k ,c k versus k. 
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3ending about 
the major axis 
of cross section 
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for the squares and products of 
v and w« 

Figure 6.- w Q /"b and-^p^ versus k 2 . 
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Figs. 7,12,13,14,15 
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Figure 12.- Section 
of cir- 
cular cylinder. 




Figure 13.- Section 
of ellip- 
tic cylinder. 
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Figure 15.- l/C k versus bending 
about minor axis.) 
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Figure 8.- Loading apparatus. 




Figure 9,- Cylinder 1 after failing 

in bending about the major 
axis of the section. 



Figure 10.- Cylinder 2 after failing 

in bending about the major 
axis of the section. 




Figure 11*- Cylinder 2 after failing 

in bending about the minor 
axis of the section. 
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Pigs. 15,17,18,19 
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Figure 18.- c versus r/ s for 

\/r=l.O according 
to tests with isotropic circular 
cylinders. 
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Figure 16.- versus n for 

different k 2 . 
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Figure 19.- c versus l/r accord- 
ing to tests on iso- 
tropic circular cylinders. 



Figure 17.- c versus r/ s . accord- 
ing to tests with 
isotropic circular cylinders. 



